Signatures of Noise Enhanced Stability in Metastable States 
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The lifetime of a metastable state in the transient dynamics of an overdamped Brownian particle 
is analyzed, both in terms of the mean first passage time and by means of the mean growth rate 
coefficient. Both quantities feature non monotonic behaviors as a function of the noise intensity, and 
are independent signatures of the noise enhanced stability effect. They can therefore be alternatively 
used to evaluate and estimate the presence of this phenomenon, which characterizes metastability 
in nonlinear physical systems. 
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Metastability is a generic feature of many nonlinear 
systems, and the problem of the lifetime of a metastable 
state involves fundamental aspects of non-equilibrium 
statistical mechanics. Metastable states, indeed, have 
been proved to play a crucial role e.g. in protein folding 
dynamics, Ising spin glasses, complex dynamics of large 
molecules at surfaces, enhancement of cellular memory, 
and in dynamics of cellular reactive oxygen species 
The problem of the lifetime of a metastable state has 
been addressed in a variety of areas, including first-order 
phase transitions, Josephson junctions, field theory and 
chemical kinetics . 

Recently, the investigation of nonlinear dynamics in 
the presence of external noisy sources led to the discovery 
of some resonance-like phenomena, among which we re- 
call stochastic resonance [j| , resonant activation Q , and 
noise enhanced stability (NES) 0, 0, ■ All these phe- 
nomena are characterized by a nonmonotonic behavior 
of some quantity as a function of the forcing noise inten- 
sity or the driving frequency, that reflects a construc- 
tive and counterintuitive effect of the noise acting on 
the nonlinear system. In particular, several theoretical 
studies have shown that the average escape time from a 
metastable state in fluctuating and static potentials has a 
non monotonic behavior as a function of the noise inten- 
sity 01 01 El- This resonance-like behavior, which contra- 
dicts the monotonic behavior predicted by the Kramers 
theory 0,0], is called the NES phenomenon: the sta- 
bility of metastable or unstable states can be enhanced 
by the noise and the average life time of the metastable 
state is larger than the deterministic decay time. Fur- 
thermore, if more realistic noise sources (such as colored 
noise with a finite correlation time) are considered, the 
value of the noise intensity at which the maximum of the 
average escape time occurs is even larger than that corre- 
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sponding to the white noise case, meaning that the NES 
effect could be easy to measure experimentally, because 
of the finite time correlations involved in any realistic 
noisy sources @. 

When considering a Brownian particle in the presence 
of a metastable potential, the NES effect is always ob- 
tained regardless on the unstable initial position of 
the particle. More precisely, two different dynamical 
regimes occur: one is characterized by a non monotonic 
behavior of the average escape time, as a function of noise 
intensity, and the other features a divergence of the mean 
escape time when the noise intensity tends to zero, im- 
plying that the Brownian particle remains trapped within 
the metastable state in the limit of small noise intensi- 
ties. The description of the transition from one dynami- 
cal regime to the other is yet an open question. 

In this Letter we analyze in more detail the divergent 
dynamical regime, and suggest an approach for detecting 
the stability of metastable states, that is alternative to 
the mean-first-passage-timc (MFPT) technique. In par- 
ticular, we find a non monotonic behavior of the MFPT 
with a minimum and a maximum as a function of the 
noise intensity for initial positions of the Brownian parti- 
cle close to the point x c where the potential shape inter- 
sects the a;— axis (see the inset of Fig.0. In this regime 
the standard deviation of the escape time has a divergent 
behavior for small noise intensity. This means that the 
average escape time, that is a quantitative measure of 
the average life time of the metastable state, has in fact 
some statistical limitations to fully describe the stability 
of metastable states. 

To complement the analysis of the transient dynamics 
of metastable states, we then introduce a different ap- 
proach, based on the evaluation of the mean growth rate 
coefficient A as a function of the noise intensity. The A 
coefficient is evaluated by the use of a procedure similar 
to that for the calculation of the Lyapunov exponent in 
stochastic systems, e.g. we consider the evolution of the 
separation 5x(t) between two neighboring trajectories of 
the Brownian particle d El 13 
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also the A coefficient displays a non monotonic behavior 
(with a clear minimum) , as a function of the noise inten- 
sity, thus representing an independent way for detecting 
and estimating the NES effect. 

The starting point of our analysis is a Brownian par- 
ticle in one spatial dimension, obeying the following 
Langevin equation 

* = -f^M + ^(*), (i) 

where £(£) is a white Gaussian noise source with zero 
mean and (5-correlated in time ( (£(t)) = and (£(£)£(£ + 
t)) = S(t)), and U(x) = 0.3x 2 — 0.2x 3 is a cubic potential, 
whose shape is shown in the inset of Fig.^ The potential 
profile has a local stable state at x s = 0, an unstable state 
at x u = 1, and intersects the x-axis at x c = 1.5. 

After fixing a given target position xp > x c , the MFPT 
t(xo,xf) (the average time for a particle starting from 
an initial position xo to reach xf) is given by the closed 
analytical expression 0, 0] 

t{xo,x f ) = 1 r e 2u ^ f e~ 2u ^ dydz , (2) 

^ J Xq J — oo 

where u(x) is the dimensionless potential profile, ob- 
tained by normalizing U(x) to the noise intensity D 
(u(x) = ^jf^). The double integral J5J) can be evaluated 
(in part analytically, in part numerically), giving rise to 

T (x ,x F ) = ^ f F e 2u ^G(z)dz, (3) 

where G{z) = 0.6046 e - z [I_ 1/3 (z) + I 1/3 (z)] - 
\ 2 F 2 (I,l;|,|;-2z) + / 2 e- 2 "(«) dy, z = 1/(101?), 
I n (z) is the modified Bessel function of the first kind 
and p Fq{a\, a?; b\, b%\ z) is the generalized hypcrgeomct- 
ric function. 

In Fig. Q] the evaluation of expression for xf = 
2.2 and for different initial positions xq (as sketched in 
the legend) is reported. Two different regimes can be 
observed, depending on the initial state of the particle: 

(i) the NES effect for all initial conditions x > x c , and 

(ii) the divergent regime for x u < x < x c , where x u = 1 
is the location of the relative maximum of the potential 
profile. 

In this latter regime, for all initial positions xo smaller 
than (but sufficiently close to) x c , the MFPT displays a 
non monotonic behavior with a minimum and a maxi- 
mum. For very low noise intensities, the Brownian parti- 
cle is trapped into the potential well, as a consequence of 
the divergence of the MFPT in the limit D —> 0. For in- 
creasing noise intensity, the particle can escape out more 
easily, and the MFPT decreases. As the noise intensity 
reaches a value D m AU = 0.1 (corresponding to the po- 
tential barrier height), the concavity of the MFPT curves 



changes. Close to such a noise intensity, the escape pro- 
cess of the Brownian particle is slowed down, due to the 
fact that the probability to reenter the well is increased. 
At higher noise intensities one recovers a monotonic de- 
creasing behavior of the MFPT. In summary, from in- 
spection of Fig. ^ the behavior of MFPT vs. D goes 
with continuity from a monotonic divergent behavior to a 
non monotonic finite behavior (typical NES effect), pass- 
ing through a non monotonic divergent behavior with a 
minimum and a maximum. 
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FIG. 1: MFPT t(xo,xf) evaluated from expression vs. 
the noise intensity D, for xf = 2.2 and for different initial 
positions xo (see the legend). The NES effect is observed for 
xq > x c , while the divergence regime corresponds to for x u < 
xo < x c . Inset: the cubic potential U(x) with the metastable 
state at x = 0. The arrow indicates the intersection point 
x c — 1.5 between the potential curve and the horizontal axis. 



In the following we will focus on this last dynamical 
regime. First of all, we compare the theoretical results 
of Fig. ^with direct numerical simulations of Eq. QJ. 
Namely, we numerically integrate Eq. Q with different 
initial conditions xq x c . For each initial condition (and 
for each value of the noise intensity D) , the integration 
is performed over an ensemble of Nr = 350, 000 different 
realizations of the white noise process The target 

position xp — 2.2 is selected, and the MFPT (calculated 
as the ensemble average of the first passage times through 
xp for different noise realizations) is reported in Fig. [21 
together with the corresponding theoretical curves. The 
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FIG. 2: Numerical estimation of the MFPT (symbols) vs. 
D for xo — 1.46 (filled circles), xq = 1.47 (empty circles), 
and Xo = 1.48 (filled triangles). The solid curves are the 
corresponding estimations of the MFPT from Eq. Inset: 
standard deviation a of MFPT's corresponding to different 
noise realizations vs. D for the same initial positions xq. 



agreement between theoretical predictions and numerical 
simulations is very good. For low values of D, however, 
we cannot reproduce the theoretical curves because of the 
finiteness of the integration time (T max = 20,000 a.u.) 
and of the number of realization. The MFPT evaluated 
in this range of noise intensities, therefore, tends to the 
deterministic escape time. 

The numerical study of Eq. (JJJ allows also to calcu- 
late the standard deviation a of the set of first passage 
times obtained for different noise realizations. The re- 
sults are shown in the inset of Fig. [5] for the same values 
of Xq, where a divergent behavior of cr(D) is visible in the 
limit D — > 0. Such a feature confirms that the only in- 
formation on the MFPT is not sufficient to fully unravel 
the statistical properties of this dynamical regime, and 
motivated our search for a complementary approach. 

This is done by monitoring the properties of the mean 
growth rate coefficient A. Let 8x$ = Sx(t = 0) < 1 be the 
initial separation of two neighboring Brownian particles 
subjected to the same noise process £(t). By linearization 
of Eq. the evolution of the particle separation Sx(t) 
is given by 



5x(t) = - 



d 2 U(x) 
dx 2 



5x(t) = \ l (x,t)Sx(t) , 



(4) 



and allows for the definition of an instantaneous growth 
rate Xi(x,t). It is important to stress that, in Eq. 

d Jgp is calculated onto (and as so, it is a function of) the 
noisy trajectory x[£(£)] The growth rate coefficient 
Aj (for the ith noise realization), is then defined as the 
long-time average of the instantaneous A^ coefficient over 



t(x ,x f ) [F 



A; 



t(x ,x f ) 



t(x ,xf) 



Xi(x, s)ds. 



(5) 



Notice that, in the limit t(xq,xp) — > oo, Eq. (JSJ) coin- 
cides formally with the definition of the maximum Lya- 
punov exponent, and that, therefore, our Aj coefficient 
has the meaning of a finite-time Lyapunov exponent, 
since we are interesting in characterizing a transient dy- 
namics. The mean growth rate coefficient (MGRC) A is 
then obtained by ensemble averaging the A, coefficients 

over the Nr different noise realizations (A = )• 

In our simulations of Eq. Q , an integration time step 
of dt = 0.0001, and an initial condition 5xq = 0.01 are 
used. The instantaneous growth rate A, (a;, t) is calculated 
over N T small sub-intervals of the trajectory, each of them 
corresponding to a time interval r = bdt. 

Instantaneous Growth Rate Coefficient 
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FIG. 3: Xi(x,t) (see text for definition) vs. time, for xo = 1.3, 
if = 3.0, and D = 1. 

In Fig. |3 we report the behavior of Xi(x,t) as a func- 
tion of time for D = 1. The Xi(x, t) coefficient is negative 
for most of the time, due to noise-induced local stability 
of metastable state. As the Brownian particle reaches 
the point at which the potential changes its concavity, 
Xi(x,t) becomes positive. As for the A coefficient, fixing 
xp = 20.0, it displays a non monotonic behavior with a 
minimum as a function of noise intensity (see Fig. 0b) , 
though being always positive. This reflects the fact that 
the Brownian particles always escape in average from the 
metastable state, but the non monotonic behavior of the 
MGRC marks the presence of the NES effect. In particu- 
lar, the closer is the initial position of the particles to x u , 
the smaller is the MGRC. Notice that Xi(x, t) is propor- 
tional to the trajectory of the Brownian particle because 
the potential is a cubic one. Any metastable state, how- 
ever, can be described through a local cubic potential 
even if the real potential has other local or global stable 
states. 

Fig. El is a comparative plot of the numerically calcu- 
lated MFPT (a) and the corresponding MGRC (b) vs. 
D for various initial positions of the Brownian particles. 
The MFPT's plot [Fig. g| (a)] shows the divergent be- 
havior as a maximum, which is shifted towards lower 
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FIG. 4: (a) Numerical estimation of the MFPT for various 
initial positions (as in the legend) and (b) the corresponding 
MGRC vs. the noise intensity D. The average has been 
performed over 20, 000 realizations. The absorbing boundary 
is if = 20.0, and the maximum waiting time used in the 
simulations is T max = 20, 000. 



values of the noise intensity. Because of the finitcness 
of the ensemble of particles considered in our numeri- 
cal experiments, observation of the divergence in time 
is prevented in our simulations, and for very low noise 
intensities (D — > 0) the deterministic escape time is re- 
trieved. As it can be seen in Fig. 01 (b), the maximum 
of the MFPT curve is reflected by a local minimum in 
the MGRC shape, which is however slightly shifted to- 
wards lower values of the noise intensity. This is be- 
cause, at low noise intensities, not all the particles coming 
back into the potential well reach positions around the 
metastable state with positive concavity, that instead will 
be attained by particles experiencing larger noise driv- 
ings. These latter ones contribute to lower down the 
value of A, due to the negative contributions of the Aj's 
that corresponds to such positions (see Fig. |2J . 

In conclusion, we investigated the average escape time 
from a metastable state in a cubic potential profile for 
different initial unstable positions of the Brownian parti- 
cles. The two introduced measures (the behavior of the 
MFPT and of the MGRC) furnish suitable tools to de- 
tect the lifetime of metastable states in the presence of 
external noise, and, as so, they can be alternatively used 



in many relevant circumstances, such as understanding 
activation processes in complex systems characterized by 
energy landscape 0. 
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